home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dptsv.z / dptsv
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPTTTTSSSSVVVV((((3333FFFF))))                                                            DDDDPPPPTTTTSSSSVVVV((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPTSV - compute the solution to a real system of linear equations A*X =
  10.      B, where A is an N-by-N symmetric positive definite tridiagonal matrix,
  11.      and X and B are N-by-NRHS matrices
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DPTSV( N, NRHS, D, E, B, LDB, INFO )
  15.  
  16.          INTEGER       INFO, LDB, N, NRHS
  17.  
  18.          DOUBLE        PRECISION B( LDB, * ), D( * ), E( * )
  19.  
  20. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  21.      DPTSV computes the solution to a real system of linear equations A*X = B,
  22.      where A is an N-by-N symmetric positive definite tridiagonal matrix, and
  23.      X and B are N-by-NRHS matrices.
  24.  
  25.      A is factored as A = L*D*L**T, and the factored form of A is then used to
  26.      solve the system of equations.
  27.  
  28.  
  29. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  30.      N       (input) INTEGER
  31.              The order of the matrix A.  N >= 0.
  32.  
  33.      NRHS    (input) INTEGER
  34.              The number of right hand sides, i.e., the number of columns of
  35.              the matrix B.  NRHS >= 0.
  36.  
  37.      D       (input/output) DOUBLE PRECISION array, dimension (N)
  38.              On entry, the n diagonal elements of the tridiagonal matrix A.
  39.              On exit, the n diagonal elements of the diagonal matrix D from
  40.              the factorization A = L*D*L**T.
  41.  
  42.      E       (input/output) DOUBLE PRECISION array, dimension (N-1)
  43.              On entry, the (n-1) subdiagonal elements of the tridiagonal
  44.              matrix A.  On exit, the (n-1) subdiagonal elements of the unit
  45.              bidiagonal factor L from the L*D*L**T factorization of A.  (E can
  46.              also be regarded as the superdiagonal of the unit bidiagonal
  47.              factor U from the U**T*D*U factorization of A.)
  48.  
  49.      B       (input/output) DOUBLE PRECISION array, dimension (LDB,N)
  50.              On entry, the N-by-NRHS right hand side matrix B.  On exit, if
  51.              INFO = 0, the N-by-NRHS solution matrix X.
  52.  
  53.      LDB     (input) INTEGER
  54.              The leading dimension of the array B.  LDB >= max(1,N).
  55.  
  56.      INFO    (output) INTEGER
  57.              = 0:  successful exit
  58.              < 0:  if INFO = -i, the i-th argument had an illegal value
  59.              > 0:  if INFO = i, the leading minor of order i is not positive
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDPPPPTTTTSSSSVVVV((((3333FFFF))))                                                            DDDDPPPPTTTTSSSSVVVV((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              definite, and the solution has not been computed.  The
  75.              factorization has not been completed unless i = N.
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.